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Abstract 

We apply recently developed effective field theory nuclear models in mean field 
approximation (parameter sets Gf and G2) to describe ground-state properties of 
nuclei from the valley of /^-stability up to the drip lines. For faster calculations of open- 
shell nuclei we employ a modified BCS approach which takes into account quasi-bound 
levels owing to their centrifugal barrier, with a constant pairing strength. We test 
this simple prescription by comparing with available Hartree-plus-Bogoliubov results. 
Using the new effective parameter sets we then compute separation energies, density 
distributions and spin-orbit potentials in isotopic (isotonic) chains of nuclei with magic 
neutron (proton) numbers. The new forces describe the experimental systematics 
similarly to conventional non-linear a — to relativistic force parameters like NL3. 



PACS: 21.60.-n, 21.10.Dr, 21.30.Fe 
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1 Introduction 



The relativistic field theory of hadrons known as quantum hadrodynamics (QHD) has be- 
come a very useful tool for describing bulk and single-particle properties of nuclear mat- 
ter and finite nuclei in the mean field approximation fl],@,|],[|]. Compared with the non- 
relativistic approach to the nuclear many-body problem, the relativistic model explicitly 
includes the mesonic degrees of freedom and treats the nucleons as Dirac particles. At the 
mean field (Hartree) level, nucleons interact in a relativistic covariant way by exchanging 
virtual mesons: an isoscalar-vector u meson, an isoscalar-scalar a meson and an isovector- 
vector p meson. With these ingredients the mean field treatment of QHD automatically 
takes into account the spin-orbit force, the finite range and the density dependence of the 
nuclear force. Adjusting some coupling constants and meson masses from the properties 
of a small number of finite nuclei, the relativistic mean field (RMF) model produces ex- 
cellent results for binding energies, root-mean-square radii, quadrupole and hexadecapole 
deformations and other properties of spherical and deformed nuclei [|5]||. 

The original linear a — u model of Walecka |7] was complemented with cubic and quartic 
non-linearities of the a meson || (non-linear a — uo model) to improve the results for the in- 
compressibility and for finite nuclei. Since these models were proposed to be renormalizable, 
the scalar self-interactions were limited to a quartic polynomial and scalar-vector or vector- 
vector interactions were not allowed. Recently, and inspired by effective field theory (EFT), 
Furnstahl, Serot and Tang [|9|JT0f| abandoned the idea of renormalizability and extended the 
RMF theory by including other non-linear scalar-vector and vector-vector self-interactions 
as well as tensor couplings HISHMIII- 

The EFT Lagrangian has an infinite number of terms since it contains all the non- 
renormalizable couplings consistent with the underlying QCD symmetries. Therefore, it is 
mandatory to develop a suitable scheme of expansion and truncation. At normal nuclear 
densities the scalar ($) and vector (W) meson fields are small compared with the nucleon 
mass (M), and they vary slowly with position in finite nuclei. This indicates that the ratios 
$/M, W/M, |V$|/M 2 and (VH^I/M 2 can be used as the expansion parameters. With the 
help of the concept of naturalness, it is then possible to compute the contributions of the 
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different terms in the expansion and to truncate the effective Lagrangian at a given level of 
accuracy H JT0| , p!2| , p!3|| . None of the couplings should be arbitrarily dropped out to the given 
order without a symmetry argument. 

References [ I0| , |I2 ,|I3| have shown that it suffices to go to fourth order in the expansion. 
At this level one recovers the standard non-linear a—uj model plus a few additional couplings, 
with thirteen free parameters in all. These parameters have been fitted (parameter sets Gl 
and G2) to reproduce some observables of magic nuclei |TtJ. The fits display naturalness (i.e., 
all coupling constants are of the order of unity when written in appropriate dimensionless 
form), and the results are not dominated by the last terms retained. This evidence confirms 
the utility of the EFT concepts and justifies the truncation of the effective Lagrangian at 
the first lower orders. 

Recent applications of the models based on EFT include studies of pion-nucleus scatter- 
ing [THl and of the nuclear spin-orbit force [ITS] , as well as calculations of asymmetric nuclear 



matter at finite temperature with the Gl and G2 sets (T^|. In a previous work |L7| we have 
analyzed the impact of each one of the new couplings introduced in the EFT models on the 
nuclear matter saturation properties and on the nuclear surface properties. In Ref. |TB[ we 
have looked for constraints on the new parameters by demanding consistency with DBHF 
calculations and the properties of finite nuclei. During the last years a large amount of 
work has been devoted to measuring masses of nuclei far from stability ||19|| . This body of 
experimental data has been used as a benchmark to test the predictions of the currently 



existent (relativistic and non-relativistic) nuclear effective forces [20]. This fact motivates 
us to investigate in the present work the behaviour of the parameter sets Gl and G2 derived 
from EFT in regions far from the stability line. 

To study ground-state properties of spherical open-shell nuclei one has to take into 
account the pairing correlations. Relativistic mean field calculations near the /3-stability 
line have usually included pairing in a constant gap BCS approximation [||,^l],[^], with the 
gaps fitted to empirical odd-even mass differences. This approach works properly when the 
main effect of the pairing correlations is a smearing of the Fermi surface. Since the BCS 
pairing energy diverges for large momenta, a cut-off has to be introduced in the pairing 
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channel to simulate phenomenologically the finite range of the particle-particle force. The 
limitations of this simple BCS method appear when one deals with nuclei far from the in- 
stability line. Close to the drip lines the Fermi level falls near the particle continuum and it 
is known that the BCS model does not provide a correct description of the coupling between 
bound and continuum states p3| , ^4| . In the non-relativistic framework this difficulty was 
overcome by the unified description of the mean field and the pairing correlations provided 
by the Hartree-Fock-Bogoliubov (HFB) theory |2^J2^], with Skyrme [^3|,^4J] or Gogny forces 



The same unified treatment was developed by Kucharek and Ring |28j in the relativistic 
framework. However, a quantitative description of the pairing correlations in nuclei cannot 
be achieved with relativistic mean field parametrizations because the meson exchange forces 
are not properly adapted to large momentum transfer p8|j29|| . Later, Ring and coworkers 
PU| , PU| , PT|P^ ] have used the RMF interaction for the particle-hole channel plus the pairing 
part of the Gogny force |27|] (with the D1S parameters f33fl ) for the particle-particle channel, 
in relativistic Hartree-plus-Bogoliubov (RHB) calculations. Other authors have employed a 
density-dependent zero-range pairing force instead of the Gogny pairing force P^JS^] . 

Recent calculations with non-relativistic Skyrme forces and a zero-range force in the 
particle-particle channel have shown that a BCS approach is able to provide a good qualita- 
tive estimate of the drip lines if some quasibound states due to their centrifugal barrier (plus 
the Coulomb barrier for protons) are included in the calculation [^,^,^]. In this work we 
will use a similar BCS approach with quasibound states to approximately take into account 
the effects of the continuum contributions near the drip lines. We will employ a constant 
pairing strength which can be considered as a simplification of the zero-range pairing force 
and which gives similar results to those obtained with a delta force for spherical nuclei |4(| . 

The paper is organized as follows. We summarize the mean field approximation to the 
EFT nuclear model in the second section. In the third section we describe our modified BCS 
approach with quasibound states, and perform some calculations to test its possibilities 
and limitations by comparing with Bogoliubov results available from the literature. The 
fourth section is devoted to the detailed study with the EFT parametrizations Gl and 
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G2 of properties such as separation energies, particle densities and spin-orbit potentials of 
nuclei belonging to chains of isotopes (isotones) with magic proton (neutron) number. Our 
conclusions are laid in the last section. 



2 Relativistic mean field approach from effective field 
theory 

The effective field theory approach to QHD has been developed in the recent years. The 
theory and the equations for nuclear matter and finite nuclei can be found in the literature 
j|,|9],[l0| and here we shall only outline the formalism. We start from Ref. where the 
field equations were derived from an energy density functional containing Dirac baryons 
and classical scalar and vector mesons. This functional can be obtained from the effective 
Lagrangian in the Hartree approximation, but it can also be considered as an expansion in 
terms of the ratios of the meson fields and their gradients to the nucleon mass of a general 
energy density functional that contains the contributions of correlations within the spirit of 
density functional theory PJTO . 



According to Refs. @,|n| the energy density for finite nuclei can be written as 
= S>a{ -ia-V + /3(M -$) + W + ±t 3 R+^^A 
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where the index a runs over all occupied states y? a (r) of the positive energy spectrum, 



$ = g s <p Q (r), W = g v Vo(r), R = g p b (r) and A = eA (r). Variation of the energy density 
([J) with respect to ip^ and the meson fields gives the Dirac equation fulfilled by the nucleons 
and the meson field equations, which are solved self-consistently by numerical iteration. We 
refer the reader to Ref. [|TI]] for the expressions of the variational equations. 

The terms with g 1) A, f3 s and (3 V take care of effects related with the electromagnetic 



structure of the pion and the nucleon (see Ref. Specifically, the constant g 7 concerns 

the coupling of the photon to the pions and the nucleons through the exchange of neutral 
vector mesons. The experimental value is g 2 /47r = 2.0. The constant A is needed to 
reproduce the magnetic moments of the nucleons. It is defined by 

A = ^A p (l + T3) + ^A„(l-r 3 ), (2) 

with A p = 1.793 and X n = —1.913 the anomalous magnetic moments of the proton and 
the neutron, respectively. The terms with f3 s and /3 V contribute to the charge radii of the 
nucleon |TjJ. 



In this work we will employ the EFT parameter sets Gl and G2 of Refs. @,|10]. The 



masses of the nucleon and the us and p mesons take their experimental values: M = 939 
MeV, m v = 782 MeV and m p = 770 MeV. The thirteen parameters m s , g s , g v , g p , rji, 
^2, f] P i «4, Co, fv, Oil and ai were fitted by a least-squares optimization procedure to 
twenty-nine observables (binding energies, charge form factors and spin-orbit splittings near 
the Fermi surface) of the nuclei 16 0, 40 Ca, 48 Ca, 88 Sr and 208 Pb, as described in Ref. |TTJ . 
The constants (3 S , /3 V and f p were then chosen to reproduce the experimental charge radii of 
the nucleon. The fits yielded two best, distinct parameter sets (Gl and G2) with essentially 
the same x 2 value [[UJ . 



We report in Table 1 the values of the parameters and the saturation properties of Gl and 
G2. One observes that the fitted parameters differ significantly between both interactions. 
For example, G2 presents a positive value of /%, as opposite to Gl and to many of the most 
successful RMF parametrizations, such as the NL3 parameter set [41]. Formally a negative 



value of K4 is not acceptable because the energy spectrum then has no lower bound J42 . 
Furthermore, the wrong sign in the $ 4 coupling constant may cause troubles in obtaining 
stable solutions in light nuclei like 12 C. We note that the value of the effective mass at 



saturation M^/M in the EFT sets (~ 0.65) is somewhat larger than the usual value in the 
RMF parameter sets (~ 0.60). This fact is related with the presence of the tensor coupling 
f v of the ui meson to the nucleon, which has an important bearing on the spin-orbit force 



One should mention that the EFT perspective also has been helpful to elucidate the 
empirical success of the usual non- linear a — uo models that incorporate less couplings (just 
up to cubic and quartic self- interactions of the scalar field): the EFT approach accounts 
for the success of these RMF models and provides an expansion scheme at the mean field 
level and for going beyond it []4|, [lT)[|l2"| . In practice it has been seen that the mean field 
phenomenology of bulk and single-particle nuclear observables does not constrain all of the 
new parameters of the EFT model unambiguosly. That is, the constants of the EFT model 
are underdetermined by the observables currently included in the fits and different parameter 
sets with low \ 2 (comparable to Gl and G2) can be found JM,|T2|Jr3,[T3 . However, the extra 
couplings could prove to be very useful for the description of further observables. Indeed, 
for densities above the normal saturation density, and owing to the additional non-linear 
couplings, the EFT models are able |TJ| to give an equation of state and nuclear matter scalar 
and vector self-energies in much better agreement with the microscopic Dirac-Brueckner- 
Hartree-Fock predictions than the standard non-linear a — u parametrizations (the latter 
completely fail in following the DBHF trends as the nuclear density grows [ |T£| , f2"2]| ) . 

The sets Gl and G2 were fitted including centre-of-mass corrections in both the binding 
energy and the charge radius. Therefore, we will utilize the same prescription of Ref. [jnj 
in our calculations with Gl and G2. Namely, a correction 

E CM = ^ MeV (3) 

to the binding energy and a correction 

" -7 1 7 fal 2 (4) 

A(2MAE CM ) 1 ' 

to the mean-square charge radius. 
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3 The pairing calculation 



It is well known that pairing correlations have to be included in any realistic calculation of 
medium and heavy nuclei. In principle the microscopic HFB theory should be used for this 
purpose. However, for pairing calculations of a broad range of nuclei not too far from the 
/9-stability line, a simpler procedure is usually considered in which a seniority potential acts 
between time-reversed orbitals. In this section we want to discuss and test a straightforward 
improvement of this simple approximation to be able to describe in addition nuclei near the 
drip lines, at least on a qualitative level. Without the complications intrinsic to a full 
Bogoliubov calculation, our faster approximation will allow us later on to perform extensive 
calculations of chains of isotopes and isotones with the relativistic parameter sets. 



The pairing correlation will be considered in the BCS approach |25],|26|. One assumes 
that the pairing interaction f pair has non-zero matrix elements only between pairs of nucleons 
invariant under time reversal: 

{a 2 a 2 \v p ^ T \aia l ) = -G , (5) 

where |a) = \nljm) and \a) = \nlj—m) (with G > and m > 0). Most often the 
BCS calculations in the RMF model have been performed using a constant gap approach 
[p,j2"I,22|. Instead, here we choose a seniority-type interaction with a constant value of G 



for pairs belonging to the active pairing shells. 

The contribution of the pairing interaction to the total energy, for each kind of nucleon 
(neutrons or protons), is 



2 

2 



£ P air = -G\ X, Ml - n a )Y" - G X K , (6) 

where n a is the occupation probability of a state with quantum numbers a = {nljm} and 
the sum is restricted to positive values of m. One has 



2 



(e a - /i) 2 + A 2 

The Lagrange multiplier [i is called the chemical potential and the gap A is defined by 

A = Gj2K(l-n a )] 1/2 . 



(7) 
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As usual the last term in Eq. (|^) will be neglected. It is not a very important contribution 
and its only effect is a renormalization of the pairing energies |]25|j26| . 



Assuming constant pairing matrix elements (|5]) in the vicinity of the Fermi level one gets 



- = 1 ( 9 ) 

2 St. v /(£ a -/i) 2 + A2 



A, (10) 



where A is the number of neutrons or protons involved in the pairing correlation. The 
solution of these two coupled equations allows one to find ii and A. Using Eqs. (|?|) and (|8]) 
the pairing energy for each kind of nucleon can be written as 

A 2 

#pair = —q- (11) 

This simple approach breaks down for nuclei far from the stability line. The reason is 
that in this case the number of neutrons (for isotopes) or protons (for isotones) increases, 
the corresponding Fermi level approaches zero and the number of available levels above it 
is clearly reduced. Moreover, in this situation the particle-hole and pair excitations reach 
the continuum. Ref. |23j showed that if one performs a BCS calculation using the same 



quasiparticle states as in a HFB calculation, then the BCS binding energies are close to the 
HFB ones but the r.m.s. radii (i.e., the single-particle wave functions) dramatically depend 
on the size of the box where the calculation is performed. This is due to the fact that there 
are neutrons (protons) that occupy continuum states for which the wave functions are not 
localized in a region, thus giving rise to an unphysical neutron (proton) gas surrounding the 
nucleus. 



Recent non-relativistic calculations near the drip lines with Skyrme forces |38| , |39| have 
shown that the above problem of the BCS approach can be corrected, in an approximate 
manner, by taking into account continuum effects by means of the so-called quasibound 
states, namely, states bound because of their own centrifugal barrier (centrifugal-plus- 
Coulomb barrier for protons). When the quasibound states are included in the BCS calcu- 
lation (from now on a qb-BCS calculation), it is necessary to prevent the unrealistic pairing 
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of highly excited states and to confine the region of influence of the pairing potential to the 
vicinity of the Fermi level. Instead of using a cutoff factor as in Ref. |J8] , in our calculations 
we will restrict the available space to one harmonic oscillator shell above and below the 
Fermi level. 

In order to check this approach we have performed with the Gl parameter set (G n = 21/ A 
MeV, see next section) calculations of the binding energy and r.m.s. radius of the 120 Sn and 
160 Sn nuclei in boxes of sizes between 15 and 25 fm fas in the non-relativistic calculations 



of Ref. [23|1). The results taking into account the quasibound levels I/19/2, 2/5/2 and U13/2 
for 120 Sn, and li\\i% and lj'13/2 for 160 Sn, are compared in Figure 1 with the output of a 
standard BCS calculation with only bound levels. It turns out that in the qb-BCS case the 
results are essentially independent of the size of the box where the calculations are carried 
out. When the quasibound levels are included the binding energies are larger than when 
only the bound levels are taken into account, due to the damping of the pairing correlation 



caused by disregarding the continuum states in the standard BCS calculation . We also 
show in Figure 1 the results of a BCS calculation using all bound and unbound levels (i.e., 
without restricting ourselves to quasibound levels) in the considered range. It is obvious 
that in this case the results are box dependent, as the binding energy and neutron r.m.s. 
radius of 160 Sn evidence. 

Another test of the qb-BCS approach concerns the asymptotic behaviour of the particle 



densities [24]. In Figure 2 we display the radial dependence of the neutron density of Sn 



(as in Ref. [23]) calculated with the Gl parameter set in boxes of radii between 15 and 25 
fm. For large enough distances the density decreases smoothly when the size of the box 
increases (except very near of the edge, where the density suddenly drops to zero because of 
the (p a = boundary condition). This means that no neutron gas surrounding the nucleus 
has appeared. In a Bogoliubov calculation the asymptotic behaviour of the particle density 
is governed by the square of the lower component of the single-quasiparticle wave function 



corresponding to the lowest quasiparticle energy |24j. This asymptotic behaviour is displayed 
by the (almost straight) dotted line in Figure 2. It can be seen that the density obtained 
with our approach decreases more slowly than the RHB density, i.e., asymptotically the qb- 
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BCS density is not able to follow the RHB behaviour. This coincides with the conclusion of 
Ref. p4] (see Figure 19 of that work) where non-relativistic HFB densities are compared for 
large distances with the densities obtained in the qb-BCS approach with a state-dependent 
pairing []3"7|] . 

Although the qb-BCS densities do not display the right asymptotic behaviour, it was 
conjectured in Ref. p4j that such an approach could allow one to compute properties of 
nuclei much closer to the drip lines than in a standard BCS calculation. Very recently, RHB 
calculations up to the drip lines of the two-neutron separation energy S 2n for nickel isotopes 



35 1 and of the charge and neutron r.m.s. radii for tin isotopes have been carried out 



using the NL-SH parameter set plus a density-dependent zero-range pairing force. We 
have repeated these calculations with our qb-BCS method for both isotopic chains (with a 
pairing interaction strength G n = 22.5/A in the case of NL-SH). 

We display the values of the S^n separation energies for the Ni chain in Figure 3a. The 
RHB calculation predicts the drip line at the isotope 100 Ni and shows shell effects at N = 28 
and 50 (and to a minor extent at N = 70). These features are well predicted by our simpler 
qb-BCS calculation. The differences between these qb-BCS and RHB results also come in 
part from the different pairing forces used in the calculations. To investigate this point we 
show in Figure 3b the neutron pairing energy obtained in our approach [Eq. fllTD ] for the 
isotopes of the Ni chain. It vanishes at N = 28, 50 and 70, in agreement with the shell 
structure shown in Figure 3a by the S^n separation energies. The largest pairing energies are 
found in the middle of two closed shells and they are enhanced by increasing N. Figure 3b 
can be compared with the RHB values displayed in Figure 2 of Ref. ||35| . The tendencies 
are the same, though the qb-BCS pairing energies are slightly larger than in the RHB 
calculation. In Figure 4 we draw our results for the radii of the Sn isotopes, and compare 
them with the RHB values. In the case of the charge radii the agreement is excellent. The 
neutron radii obtained in our method closely follow the behaviour of the RHB neutron radii 
and the kink at N = 132 is qualitatively reproduced. 

We have furthermore computed the binding energies of nuclei of the N = 20 isotonic 



chain for which RHB results exist with the NL3 parameter set [31j|. We present the extracted 
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two-proton separation energies S% p in Table 2. The agreement between the qb-BCS and 
RHB approaches again is very good. In both models the last stable nucleus is 46 Fe, as in 
experiment. Notice that in the present case the first levels with positive energy correspond 
to those of the pf shell. Due to the Coulomb barrier all these levels become quasibound 
in our approach, and it is expected that they will lie close to the canonical levels. This 
explains the goodness of the qb-BCS energies for this isotonic chain. 

^From the previous comparisons we see that the simple qb-BCS calculation is able to 
reasonably follow the main trends of the more fundamental RHB pairing calculation. One 
can also conclude that the consideration of quasibound states in the BCS approach is, 
actually, a key ingredient to eliminate the spurious nucleon gas arising near the drip lines. 



4 Results for EFT parameter sets 

We want to analyze the ability of the Gl and G2 parameter sets based on effective field 
theory ^.HJ to describe nuclear properties far from the stability line, i.e., far from the 
region where the parameters were fitted. To our knowledge such calculations have not been 
explored so far. We will contrast the results with experiment and with those predicted by 
the NL3 set, that we take as one of the best representatives of the usual RMF model with 
only scalar self-interactions. 

As indicated, we shall use a schematic pairing with a state-independent matrix element 
G T = C T /A, where C T is a constant and r = n,p for neutrons or protons, respectively. We 
fix the constant C n for neutrons by looking for the best agreement of our calculation with the 
known experimental binding energies of Ni and Sn isotopes. Similarly, we determine C p for 
protons from the experimental binding energies of the isotones of N = 28 and N = 82. The 
values obtained from this fit are C n = 21 MeV and C p = 22.5 MeV for the Gl set, C n = 19 
MeV and C p = 25 MeV for G2, and finally C n = 20.5 MeV and C p = 23 MeV for NL3. 
Figure 5 shows that the neutron and proton state-independent gaps (A n and A p ) predicted 
by our calculation with Gl are scattered around the empirical average curve 12/y/A j44j . 
A similar picture is found with the parameter sets G2 and NL3. 
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4.1 Two-particle separation energies 

In Figure 6a we present the two-neutron separation energies S2 n for the chain of Ni isotopes. 
Clear shell effects arise at N = 28 and 50. The three relativistic interactions (Gl, G2 and 
NL3) slightly overestimate the shell effect at iV = 28 as compared with the experimental 
value, which also happens in more sophisticated RHB calculations with NL3 [p0| , p2| . In 
our qb-BCS approach some disagreement with experiment is found for the N = 38 and 
iV = 40 isotopes. Again, this also occurs in the RHB calculations of Refs. with NL3. 



However, if we compare Figure 6a with the results that we have shown in Figure 3a for the 
NL-SH parameter set, we see that NL-SH achieves a better agreement with experiment for 
these N = 38 and N = 40 isotopes. 

We stop our calculation towards the neutron drip line when the two-neutron separation 
energy vanishes or when the neutron chemical potential becomes positive. The fact that S2 n 
is not always zero at the drip line is connected with the quenching of the shell structure with 



N, which is a force-dependent property [24]. This effect is illustrated in Figure 25 of Ref. 
p4| for HFB calculations with different non-relativistic forces. We find similar situations 
with the considered relativistic sets in our qb-BCS calculations of separation energies. In 
the case of the Ni isotopes we reach the drip line at N = 66 with the Gl and NL3 sets 
and at N = 68 with the G2 set. This agrees nicely with the value N = 66 obtained in 
HFB calculations with the Skyrme forces SIII @g| and SkP @. For NL-SH our qb-BCS 



scheme predicts the drip line at N = 72 (see Figure 3a), the same value found in the RHB 
calculations of Ref. . 

In Figure 7a we display our qb-BCS results for the two-neutron separation energies of 
the Sn isotopic chain. In Ref. |32| it was claimed that pure BCS calculations in the constant 
gap approach (with NL3) are not suitable for the Sn isotopes. We observe in Figure 7a that 
below N = 60, as one moves towards N = 50, some discrepancies with the experimental 
values appear, which also arise in the RHB calculations ||32|| . The three forces slightly 
overestimate the shell effect at N = 82 (as the RHB results of Refs. [0,0 for NL3). We 
have computed Sn isotopes up to A = 176, when S2 n vanishes for NL3 (in good agreement 



with RHB results for NL-SH [36] and HFB results for the Skyrme force SkP [23]). For Gl 
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and G2 we find that S^n does not yet vanish at iV = 126, and it is not possible to increase 
the neutron number due to the shell closure at N = 126 (the neutron chemical potential 
becomes positive for the N = 128 isotope). This means that the quenching of the shell 
effect at N = 126 for NL3 (and NL-SH) is larger than for the Gl and G2 parameter sets. 

Our calculated S^n energies for Pb isotopes are shown in Figure 8. The experimental 
shell effect at iV = 126 is reasonably well reproduced by Gl, G2 and NL3. The drip 
line is found at N = 184 with a non-vanishing two-neutron separation energy, as in the 
calculations performed with the extended Skyrme force SLy4 in Ref. flSSH , where a similar 
approach (quasibound states and a state-dependent gap) to ours was used. The relatively 
large shell effect found at iV = 184 means that there is no quenching for this magic number 
in our qb-BCS approximation for the studied parameter sets. Indeed, to verify this point a 
full RHB calculation should be performed. 

To analyze the proton pairing we have studied the two-proton separation energies in 
chains of isotones of N = 28 (Figure 9a) and N = 82 (Figure 10a). In the case of N = 28, 
shell gaps appear at Z = 20 and Z = 28. For Z = 20 the predicted gap is larger than in 
experiment. For Z = 28, Gl and G2 agree better than NL3 with experiment. S 2p vanishes 
at Z = 30 for NL3, whereas it vanishes at Z = 32 for Gl and G2. The isotones of iV = 82 
display a clear shell effect at Z — 50, in agreement with the non-relativistic calculation of 
Ref. p8[ . It is slightly larger for G2 than for Gl and NL3. Experimental information for 
this shell effect is not available. NL3 would predict another shell effect at Z — 58, which 
does not appear experimentally. The effect is less pronounced in Gl and it does not show up 
in G2. The three forces indicate that the proton drip line is reached after the 156 W isotope, 
in agreement with experimental information [[16] . 

Figures 11a and lib show, respectively, the calculated S 2p separation energies for the 
N = 50 and iV = 126 isotone chains. Note that we did not use any information about these 
nuclei in our fit of the G p pairing strength. For N = 50 the set G2 follows the experimental 
data very well, specially for the larger Z. The trend of Gl and NL3 is only a little worse. 
The proton drip line is located at 100 Sn in the three parametrizations, in good accordance 
with experiment. The quenching of the shell effect at Z = 50 is larger for G2 than for 
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Gl and NL3. The available data for two-proton separation energies oi N = 126 isotones 
are reasonably well estimated by the relativistic sets. However, the trend of NL3 is worse 
than that of Gl and G2. It would then be very interesting to perform RHB calculations 
of this chain to confirm the behaviour of NL3. The last nucleus of the chain stable against 
two-proton emission is 218 U according to Gl and NL3, and 220 Pu according to G2. The 
three sets predict a shell effect at Z = 92, though it is relatively quenched for G2. 



4.2 One-particle separation energies 

We have computed one-neutron (one-proton) separation energies for Ni and Sn isotopes 
(for N = 28 and N = 82 isotones). The results are displayed in Figures 6b and 7b (9b 
and 10b), respectively. To deal with odd mass number nuclei we have used a spherical 
blocking approximation. One pair of conjugate states \a) and \a) is blocked, i.e., taken 
out of the pairing scheme [ 25|pB |. In the spherical approximation one replaces the blocked 



single-particle state by an average over the degenerate states in its j shell. This way the 
rotational and time-reversal invariance of the many-body system is restored in the intrinsic 
frame [H71 . In this approach the contribution of the j shell that contains the blocked state 



to the number of active particles and the pairing energy is 

Aj = (2j- 1K + 1, (12) 

£ pairj = -G(2j-l)[n j (l-n j )) 1 / 2 , (13) 

respectively. The remaining active shells contribute in the usual manner [Eqs. @ and fllPP]. 
Due to rearrangement effects, blocking the single-particle state with smallest quasiparticle 



energy E a = y(s a — /i) 2 + A 2 in the even A — 1 nucleus, does not necessarily lead to the 
largest binding energy of the odd A nucleus. Therefore, in some cases one has to repeat 
the calculation blocking in turn the different single-particle states that lie around the Fermi 
level to find the configuration of largest binding energy ]23||j7 , 47 . 



The one-neutron (one-proton) separation energies lie over two different curves for even 
and odd neutron (proton) number. For Ni isotopes (Figure 6b) the three parameter sets 
Gl, G2 and NL3 reproduce reasonably the experimental values. The shell effect at N = 28 
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is, again, overestimated by the three forces. The heaviest Ni isotope stable against one- 
neutron emission is found at iV = 55 with NL3 and Gl and at N = 57 with G2. For 
Sn isotopes (Figure 7b) the shell effect at N = 82 is slightly overestimated by the studied 
forces. The predictions of the three parametrizations are roughly similar up to N = 110, 
where the behaviour of NL3 starts to depart from Gl and G2 due to the large quenching of 
the N = 126 shell effect shown by NL3 as compared with Gl and G2. In our calculations, 
the odd Sn isotopes become unstable against the emission of one neutron around N = 110. 
This value is larger than the values found with the non-relativistic SkP interaction in HFB 



(N = 103) or HF+BCS (N = 101) calculations [23]. The origin of this discrepancy lies in 
the fact that the shell distribution in tin isotopes around N = 126 for SkP is rather different 
from that of the relativistic sets |32 . 



The one-proton separation energies for the isotones of N = 28 (Figure 9b) show an 
overall good agreement with the experimental data. The shell effects at Z = 20 and Z = 28 
are rather well reproduced by the forces analyzed here. Gl and G2 predict the heaviest 
nucleus stable against one-proton emission to be 57 Cu, as in experiment, while it is unstable 
in the NL3 calculation. For the isotones of N = 82 (Figure 10b) the shell effect predicted 
by NL3 and Gl at Z = 50 is similar. Again, as for S^p (Figure 10a), NL3 predicts a shell 
effect at Z = 58 which is not found experimentally, whereas for Gl this effect is clearly 
smaller and it does not appear for G2. The last stable nucleus against one-proton emission 
is 151 Tm according to the three parameter sets. 

4.3 One-body densities and potentials 

The nuclear densities for chains of isotopes of light and medium size nuclei have recently 
been studied in the RHB approximation |]3T),[31, 3~5l , ]36|| . As N grows the neutron and mass 



densities extend outwards and the r.m.s. radii and the surface thickness increase. Special 
attention has been payed to the isospin dependence of the spin-orbit interaction. The 
magnitude of the spin-orbit potential is reduced when one approaches the neutron drip 
line and, as a consequence, there is a reduction of the energy splittings between spin-orbit 
partner levels |30|j3T|j36[| . To our knowledge, for isotones such an study has only been carried 
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out in the N = 20 chain fl5I| . It is to be remarked that the EFT parametrizations Gl and 
G2 contain a tensor coupling of the u meson to the nucleon which plays a very important 
role in the spin-orbit force because there exists a trade-off between the size of this coupling 
and the size of the scalar field P~5] , [TT| ] . 

In Figures 12a and 12b we display, respectively, the neutron and proton densities of some 
N = 28 isotones from Z = 16 to Z = 32 as predicted by the G2 set in our qb-BCS approach. 
Figures 13a and 13b show the results for some N = 82 isotones from Z = 40 to Z = 70. 
Since N is fixed in each isotonic chain, the spatial extension of the neutron densities is very 
similar for the different nuclei of the chain. In any case, as one goes from the lightest to 
the heaviest isotone of the chain, the neutron densities tend to be depressed in the interior 
region and their surface thickness (90%-10% fall-off distance) shows a decreasing tendency. 
The proton densities of the isotones exhibit a strong dependence on Z: by adding more 
protons they are raised at the interior and their surface is pushed outwards. For iV = 28 
the surface thickness of p p remains roughly constant up to Z = 28 and increases for heavier 
isotones as a consequence of the growing occupation of the I/7/2 shell. At the origin the 
proton densities show a bump when Z > 20 because the 2s\/2 level is occupied. The Z = 16 
isotone shows a dip at the centre due to, precisely, the emptiness of this 2s\ji level. For the 
considered nuclei of N = 82, the proton densities have an approximately constant surface 
thickness and present a hole at the centre owing to the Coulomb repulsion. In Figure 14 we 
display the neutron and proton r.m.s. radii of the N = 82 isotonic chain obtained with G2 
and NL3. It turns out that the predictions of both sets are very similar. The proton radii 
increase uniformly with Z, similarly to the behaviour found for N = 20 isotones with NL3 
and the RHB scheme in Ref. |3T|. The neutron radii remain roughly constant with Z. They 
just show a slight decrease with increasing Z till Z ~ 50 and slightly increase afterwards. 
This behaviour may be related with the shell effect for protons at Z — 50. 

The spin-orbit interaction is automatically included in the RMF approximation. It 
appears explicitly when the lower spinor of the relativistic wave function is eliminated in 
favour of the upper spinor. This way one obtains a Schrodinger-like equation with a term 
Vso(r) that has the structure of the single-particle spin-orbit potential. Including the 
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contribution of the tensor coupling of the to meson, the spin-orbit term reads |T5|,|30|| 



Hso = ^V so (r)L-S, (14) 
T , . . M 2 l/d$ dW\ n M\dW 

V so (r) = =2- — + — )+2f v = - — , 15 
M r \ar ar J M r ar 

where M = M — |($ + W). We have checked numerically that the contribution to the 
spin-orbit potential of the f p tensor coupling of the p meson is very small, even when one 
approaches the drip lines. Hence we have not written this contribution in Eq. (|l~5l). 

The spin-orbit potential (|T5| ) for some lead isotopes computed with G2 and NL3 is 
displayed in Figures 15a and 15b, respectively. As a general trend, for both G2 and NL3, 
when the number of neutrons is increased the depth of the spin-orbit potential decreases 
gradually and the position of the bottom of the well is shifted outwards, which implies a 
significant weakening of the spin-orbit interaction. The same effect arises in other isotopic 



chains in RHB calculations [|30| , |35| , p6| . Comparing the spin-orbit potentials obtained with 
the G2 and NL3 sets, one sees that they have a similar strength for all the isotopes analyzed 
and that the minima of the wells are located at similar positions (slightly shifted to larger 
values of r in G2). The higher effective mass of G2 at saturation (M^/M = 0.664) with 
respect to NL3 (M^/M = 0.595) is compensated by the tensor coupling included in G2 
(fv — 0.692). To ascertain the relative importance of the tensor coupling we have drawn 
in the insert of Figure 15a, for 228 Pb, the full potential (|i~5|) and the contribution resulting 
from setting f v = in Eq. (p~5|) . We see that the full Vsoij) is much deeper and wider. The 
maximum depth of Vso{ r ) changes from —68 MeV fm~ 2 (right scale of the insert) to —44 
MeV fm -2 when f v = 0. That is, the tensor coupling accounts for roughly one third of the 
total spin-orbit strength in the G2 parameter set. 

One expects that the weakening of the spin-orbit potential in going to the neutron drip 
line will bring about a reduction of the spin-orbit splittings 

As = £ n /,i=i-l/2 — £ nl,j=l+l/2 (16) 

of the neutron levels P0| . Figure 16 displays the energy splittings of some spin-orbit partner 
levels of neutrons for lead isotopes, obtained with the G2 and NL3 parameter sets. The 

18 



splittings predicted by G2 and NL3 are very close as a consequence of the similarity of the 
corresponding spin-orbit potentials. Partner levels with high angular momentum undergo 
some reduction in the splitting along the Pb isotopic chain, but partners with small angular 
momentum show an almost constant splitting. By comparison of their RHB results for Ni 
and Sn, the authors of Ref. |K| pointed out that the weakening of the spin-orbit interaction 
should be less important for heavier isotopic chains. Our calculations for Pb would confirm 
this statement. All the single-particle levels involved in Figure 16 are bound. Of course, 
one should not expect the results for Ae to be so reliable in our qb-BCS approach if one, 
or both, of the partner levels lies at positive energy. The reason is that the single-particle 
energies of the quasibound levels do not exactly reproduce the energies of the corresponding 
canonical states of a RHB calculation. 

In Figures 17a and 17b we show the spin-orbit potential for isotones of N = 82 from 
Z = 40 to Z = 70, for the G2 and NL3 parametrizations. Similarly to what is found for 
isotopes, the results obtained from G2 and from NL3 are comparable and the spin-orbit 
potential well Vso{ r ) moves outwards with the addition of protons, following the tendency 
of the proton density. However, for isotones we find that the behaviour of the depth of the 
spin-orbit potential well is not so monotonous: it increases when one goes from the neutron 
drip line up to the /3-stable region, while it decreases afterwards as more protons are added. 



5 Summary and conclusion 

We have analyzed the pairing properties of some chains of isotopes and isotones with magic 
Z and N numbers in the relativistic mean field approach. The study has been performed for 
the Gl and G2 parametrizations that were obtained in Ref. [TO from the modern effective 
field theory approach to relativistic nuclear phenomenology. We have compared the results 
with those obtained with the NL3 parameter set which is considered to be very successful 
for dealing with nuclei beyond the stability line. 

For accurate calculations of pairing far from the valley of /3-stability in the relativistic 
models, the relativistic Hartree-Bogoliubov approach should be applied. However, we have 
presented a simpler modified BCS approach which allows one to obtain pairing properties 
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near the drip lines fast and confidently. The method has been used previously in non- 
relativistic calculations with Skyrme forces |38| , |39|| . The key ingredient is to take into account 
the continuum contributions through quasibound levels due to their centrifugal barrier. 
To further simplify the calculations we have assumed pairing matrix elements of the type 
G = C/A instead of, e.g., a state-dependent pairing with a zero-range force. 

The considered quasibound levels are mainly localized in the classically allowed region 
and decrease exponentially outside it. This eliminates the unphysical nucleon gas which, 
near the drip lines, surrounds the nucleus when all available positive energy levels are 
included in the usual BCS approach. Normally, the quasibound levels have high angular 
momentum and lie close in energy to the corresponding RHB canonical levels. One of the 
limitations of the qb-BCS approach employed here is the fact that the nuclear density does 
not follow the asymptotic fall-off of the densities computed with the relativistic Hartree- 
plus-Bogoliubov theory. In spite of this shortcoming, we have shown by comparison with 
available RHB results that the qb-BCS approach is able to predict the position of the drip 
lines, or the behaviour of the neutron and charge radii for nuclei far from the stability line, in 
a reasonable way. Also, the obtained pairing gaps are nicely scattered around the empirical 
average 12/ \/~A. 

We have applied the qb-BCS approach to the Ni, Sn and Pb isotopic chains, and to the 
N = 28, 50, 82 and 126 isotonic chains. The two-neutron (two-proton) and one-neutron 
(one-proton) separation energies, as well as the resulting shell gaps, are similar for the three 
studied relativistic parametrizations (Gl, G2 and NL3) and in general they reproduce the 
available experimental data, at least qualitatively. The neutron and proton drip lines are 
usually reached at the same place with the three forces, though one may find a shift of one 
or two units of A among them. We have paid some attention to the quenching of the shell 
structure near the drip lines. For example, the quenching of the shell effect at iV = 126 
for Sn isotopes is larger in NL3 than in Gl and G2, while for Pb isotopes none of the 
three sets exhibits a quenching of the shell effect at N = 184 in our qb-BCS calculation. 
The EFT parametrizations Gl and G2 contain tensor couplings that are not present in the 
RMF parametrizations like NL3 and have a larger effective mass at saturation. However, 
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the predicted spin-orbit potentials along the isotopic and isotonic chains do not differ much 
from those obtained with NL3. 

Our analysis shows that the parameter sets based on EFT are able to describe nuclei far 
from the /3-stability line, after adding a phenomenological pairing residual interaction. Only 
experimental information about some magic nuclei was utilized in the fit of the constants 
of the Gl and G2 sets and, thus, the results for nuclei near the drip lines are veritable 
predictions of the model. In spite of the fact that the EFT sets include more couplings 
and parameters than the conventional RMF sets like NL3, both models reproduce the ex- 
perimental systematics with a similar quality. In fact, the studied properties away from 
the valley of ^-stability do not seem to provide further constraints on the EFT parameters, 
not even in the isovector sector. In conclusion, extended sets like Gl and G2 will serve 
almost the same purposes for normal systems as the conventional parameter sets. However, 
some of the extra parameters of the general EFT functional may be used to better describe 



regions of the equation of state at higher density or temperature fllgjlql without spoiling 
the systematics for finite nuclei. 
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Figure captions 

Figure 1. Dependence of the binding energy (left) and neutron r.m.s. radius (right) of the 
nuclei 120 Sn and 160 Sn on the size of the box used in the calculations. Solid, dashed and 
dotted lines correspond to a BCS approach including quasi-bound levels, only bound 
levels and all available levels, respectively. The results are for the Gl parameter set. 

Figure 2. Neutron density of 150 Sn for different sizes of the box used in the qb-BCS calcu- 
lations (for the set Gl). The dotted line denotes the asymptotic behaviour expected 



from a Bogoliubov calculation |24 . 



Figure 3. Part (a): the two-neutron separation energies for Ni isotopes calculated in the 



qb-BCS approach are compared with the RHB results of Ref. |35[ and with experiment. 
Part (b): the neutron pairing energy obtained in the qb-BCS approach. The results 
are for the set NL-SH. 

Figure 4. Charge and neutron r.m.s. radii of Sn isotopes in qb-BCS and RHB [[36|] calcu- 
lations performed with the NL-SH set. 

Figure 5. The state-independent pairing gaps predicted by our qb-BCS approach for Ni, 
Sn and Pb isotopes (top) and for N = 28, 50, 82 and 126 isotones (bottom). The Gl 



set has been used. The empirical average curve 12/ y A [P| is depicted by a solid line. 



Figure 6. Two-neutron (a) and one-neutron (b) separation energies for Ni isotopes com- 
puted with the qb-BCS approach for the parametrizations Gl, G2 and NL3, in com- 
parison with the experimental data. 

Figure 7. Same as Figure 6 for Sn isotopes. 

Figure 8. Two-neutron separation energies for Pb isotopes computed with the qb-BCS 
approach, in comparison with experiment. 

Figure 9. Two-proton (a) and one-proton (b) separation energies for N = 28 isotones 
computed with the qb-BCS approach, in comparison with experiment. 
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Figure 10. Same as Figure 9 for N = 82 isotones. 

Figure 11. Two-proton separation energies for N = 50 (a) and N = 126 (b) isotones 
computed with the qb-BCS approach, in comparison with experiment. 

Figure 12. Radial dependence of the neutron (a) and proton (b) densities of some N = 28 
isotones obtained with the G2 set. 

Figure 13. Same as Figure 12 for some iV = 82 isotones. 

Figure 14. Neutron and proton r.m.s. radii of N = 82 isotones obtained with the G2 and 
NL3 sets. 

Figure 15. Spin-orbit potential for some Pb isotopes obtained with the G2 set (a) and 
with the NL3 set (b). 

Figure 16. Energy splitting of some spin-orbit partner levels of neutrons in Pb isotopes, 
calculated in the qb-BCS approach for the G2 and NL3 sets. 

Figure 17. Spin-orbit potential for some N = 82 isotones obtained with the G2 set (a) 
and with the NL3 set (b). 
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Table 1: Dimensionless parameters and saturation properties of the sets Gl and G2 based 
on EFT and of the RMF set NL3. 
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Table 2: RHB and qb-BCS two-proton separation energies (in MeV) of some N = 20 isotones 
calculated with the NL3 parameter set. 
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